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ABSTRACT 

We compute the exponential disc scalelength for 686 disc galaxies with spectroscopic 
rodshifts out to redshift 5.8 based on Hubble Space Telescope archival data. We com- 
pare the results with our previous measurements based on 30000 nearby galaxies from 
the Sloan Digital Sky Survey. Our results confirm the presence of a dominating expo- 
nential component in galaxies out to this redshift. At the highest redshifts, the disc 
scalelength for the brightest galaxies with absolute magnitude between —24 and —22 
is up to a factor 8 smaller compared to that in the local Universe. This observed 
scalelength decrease is significantly greater than the value predicted by a cosmological 
picture in which baryonic disc scalelength scales with the virial radius of the dark 
matter halo. 
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low redshift ~ 1.5 suggest no disc size evolution out to these 
redshifts (e.g., see Simard et al. 1999; Ravindranath et al. 
2004; Barden et al. 2005), whereas some studies have ap- 
plied direct size measurement methods out to redshift 5 and 
beyond, and found support for disc evolution in agreement 
with the H^'^^'^ model (Bouwens et al. 2004; Oesch et al. 
2010). Ferguson et al. (2004) and Hathi et al. (2008), on 
the other hand, argued for disc size evolution following the 
model, which implies a constant galactic rotation curve 
in a cosmological picture in which baryonic disc scalelength 
scales with the virial radius of the dark matter halo (e.g., 
Mao et al. 1998). While these studies have investigated high 
redshift surveys carried out with the Hubble Space Telescope 
(HST) . their results on disc size evolution arc based on red- 
shifts based on photometric measurements or drop-out tech- 
niques combined with direct size derivation methods. Direct 
size measurements are strictly limited by the dotoctability 
of the large lower surface brightness galaxies in any given 
magnitude in the sense that they arc not able to pick out 
these gala^xies and their correct sizes in the HST images. 

To date, high redshift studies of scalelengths have found 
a factor two decrease at the mean redshift ~ 2.5 (Elmogrcen 
et al. 2005a) and a more detailed test bed for cosmological 
simulations of galaxy formation is still missing. Here, we 
present scalelengths for 686 galaxies out to redshift 5.8 based 



1 INTRODUCTION 

A versatile observable used to establish the formation his- 
tory of a galaxy disc is the exponential scalelength (Free- 
man 1970; Lin & Pringle 1987). It relates directly to the 
disc mass and dark matter angular momentum, and thereby 
it can be used to derive the physical conditions of the grav- 
itational collapse phase and galaxy merger remnants (e.g., 
Dalcanton et al. 1997; Button 2009). It also dictates the evo- 
lution of substructures that may change the morphology of 
their host discs, caused by gravitational instabilities such as 
bulges, bars, and spiral arms (Bournaud ct al. 2007). Disc 
scalelengths can also be converted into disc sizes and masses 
with the advantage that scalelength derivation is relatively 
independent on observed luminosity since it does not only 
rely on the dctcctability limit of a given surface brightness 
contour, e.g., effective radius of a galaxy, but it also relies 
on the brighter surface brightness contours on the face of 
the disc. 

Galaxy sizes scale with redshift z following three pos- 
sible scaling relations: at fixed halo circular veloc- 
ity, H^^^'^ at fixed halo mass, and no evolution, where 
H = H{z) is the Hubble parameter which scales approxi- 
mately as (1 -|- z)^^^ (Carroll et al. 1992). Observations be- 
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Figure 1. Example light profiles from the HST sample with the 
best fit exponential profile (linear fit in log-scale). Black points 
represent the radii that were considered in each fit, and for clarity, 
the error bars are not displayed for all other points. Similar plots 
for the SDSS sample can be found in Fathi et al. (2010). The 
j/-axis is given in arbitrary units. 



on archival data from major high redshift surveys carried out 
with the HST. 



2 DATA ANALYSIS 

We have analyzed all publicly available archival imaging 
data from the three deepest ever Hubble Space Telescope 
- Advanced Camera for Surveys (HST-ACS) imaging sur- 
veys, namely the Great Observatories Origins Deep Sur- 
veys (GOODS) South and North (Giavalisco et al. 2004) 
and the Hubble Ultra Deep Field (HUDF) (Beckwith et al. 
2006). Our aim is to compare scalelengths for high redshift 
galaxies observed with the HST with our previous results for 
~ 30000 r-band selected Sloan Digital Sky Survey (SDSS; 
York et al. 2000) galaxies out to redshift 0.3 (Fathi et al. 
2010; Fathi 2010). The disc structure in the SDSS sample 
was quantified using our own customised one-dimensional 
light profile fitting, that was extensively tested and shown 
to deliver disc scalelengths very similar to previously pub- 
lished scalelengths derived by two-dimensional image fitting 
and multiple component fitting algorithms (further details 
can be found in Fathi et al. 2010 and Fathi 2010). 

Images are available in the ACS F435W, F606W, 
F775W and F850LP fihers for all three high redshift HST 
surveys. We used the best quality images in the F850LP 
filter (z— band) with Signal-to-Noise ratio > 15 and select 
all galaxies (1642 in total) for which spectroscopic redshifts 
available in the literature (Barger et al. 1999; Cohen et al. 
2000; Cowie et al. 2004; Vanzella et al. 2004, 2006, 2008). 
After selecting galaxies for which reliable images and solid 
or reliable redshift measurements (denoted by quality flags 
A or B in the online catalogues) are available, we inves- 
tigated the presence of a dominant disc-like component in 
all galaxies. In the absence of a morphological classification 
catalogue for these samples, we used the diagnostic Sersic 
profile fitting to the luminosity profile of each galaxy and 
derived a mean Sersic index 1.3±0.4, in full agreement with 
high redshift observations (Toft et al. 2005; van der Wei et 
al. 2011) and cosmological galaxy formation models (Mo et 
al. 1998; Dutton 2009). The average Sersic index of 1.3±0.4 



strongly suggests the presence of a dominant exponential 
luminosity profile in our sample galaxies, often associated 
with a galactic disc. At all redshifts, approximately 25% 
have Sersic indices above 2, and less than 10% above 4, thus, 
the contamination by E/SO galaxies in our sample does not 
have a significant effect on the derived scalelengths. 

We then quantified the structure of these exponential 
discs using the routines extensively tested in Fathi et al. 
(2010) and Fathi (2010). A brief summary follows: 

First an image tile is created centred on each galaxy, fol- 
lowed by the determination of the galaxy's isophotal major- 
and minor-axis and position angle, using SExtractor (Berlin 
& Arnouts 1996). Stars and nuisance sources are masked 
out, and asymmetry measure and inclination are calculated. 
These can help identify and discard unreliable results and 
create a sample of high redshift galaxies with identical char- 
acteristics directly comparable to the SDSS. Each galaxy 
is divided into concentric rings with constant axis ratio as 
found in the previous step. A galactocentric light profile is 
created with the mean luminosity of each ring at every ra- 
dius, fitted by 7(r) = Joe'-'"^'"'*-', where /(r) is the measured 
intensity at galactocentric each radius r, 7o its value at the 
centre of the galaxy and rd the characteristic disc scale- 
length. The inner 50% of each galaxy and the regions out- 
side two galactic radii, are omitted from the fits, hence our 
fits are not affected by the presence of central bulges, cir- 
cumnuclear starbursts, noisy or sky-dominated outskirts of 
galaxies. For all galaxies the choice of the outer fitting ra- 
dius is justified by the small relative standard deviation for 
the derived average magnitudes at each radius, i.e., always 
less than one magnitude (Fathi et al. 2010). Example fits 
are shown in Fig. 1. 

We discard galaxies for which i) the inclination is 
greater than 60 degrees to avoid discs that have been di- 
luted due to projection (Fathi et al. 2010); ii) asymmetry 
measure is greater than 0.35, to remove disturbed and unre- 
laxed systems (Shade et al. 1995); in) the disc light profile 
is represented by ^ 5 data points; iv) the least square fit has 
a relative error greater than unity. Applying these criteria 
to the initial sample of 1642 galaxies with spectroscopic red- 
shifts, we isolated a sub-sample of 686 galaxies with these 
properties identical to the 30000 SDSS galaxies studied in 
(Fathi et al. 2010). The final HST sample spans out to red- 
shift 5.8, corresponding to ~ 12.5 Gyr look-back time in 
a standard ACDM cosmology where Ho = 70 km/s/Mpc, 
Q.A = 0.7, and Q.m = 0.3 (see Fig. 2). 

We combine the scalelengths derived from the SDSS 
sample below redshift 0.3 with the HST sample, and divide 
the fuU set in three redshift bins (0-0.3; 0.3-2.0; 2.0-5.8). 
The nearby redshift bin includes 30210 galaxies with abso- 
lute magnitudes in the range (—24,-14). The intermediate 
bin includes 450 galaxies with absolute magnitudes between 
—24 and —18, and the highest redshift bin includes 26 galax- 
ies all brighter than absolute magnitudes —22. We further 
divide each bin into absolute magnitude subsets of two or- 
ders of magnitude, i.e., roughly corresponding to one order 
of magnitude in total galaxy mass. Each subset is thus pop- 
ulated by at least 26 galaxies. For each redshift and abso- 
lute magnitude bin, the mean galaxy disc scalelength values 
are quantified fitting a Gaussian function to the V/Vmax 
volume-corrected distribution histogram of each bin (e.g., 
van der Kruit 1987). Figure 3 illustrates three such his- 
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Figure 2. Exponential disc scalelengths versus spectroscopic redshift over the last 12.5 Gyr. The arc of data points at the left shows 
the 30000 SDSS galaxies while the filled circles are the 686 HST galaxies at redshift < 5.8. The error bars illustrates the approximate 
scalelength uncertainties (35%) due to observational effects. The dotted curve that illustrates the upper envelope for the decrease 
predicted for the bright nearby galaxies with large scalelengths. Out to redshift ^ 1—2, the scalelength decrease is in accordance with the 
theoretical predictions, whereas for the galaxies with absolute magnitudes between —24 and —22 and at redshifts > 2, scalelengths decrease 
significantly greater than the prediction by the theoretical H~^/^ model. The trend remains unchanged when using the photometric 
redshift sample with 1472 galaxies. The SDSS sample selection following fixed stamp size around each galaxy is illustrated by the solid 
curve at the top left of this diagram (Fathi et al. 2010). 



tograms with corresponding fits for the absolute magnitude 
bin (-24,-22). 

We find that while brighter discs are typically more 
extended at any given redshift, the disc scalelength is de- 
creasing at high redshift for any given intrinsic luminosity 
bin (see Fig. 2). At redshift < 0.3, we derive mean scale- 
length of 4.8 ± 1.5 kpc for galaxies with observed absolute 
magnitudes between —24 and —22. Their counterparts in 
the redshift range between 0.3 and 2.0 have a mean scale- 
length of 1.9 ±0.6 kpc, and at redshifts between 2.0 and 5.8, 
0.6 ±0.2 kpc (see Fig. 3). Similar comparison for fainter in- 
trinsic luminosity bins does not yield equally robust results 
due to the absence of spectroscopically confirmed galaxies 
with absolute magnitude fainter than —22 in the redshift 
bin 2.0-5.8. To examine this effect for fainter galaxies, we 
make use of all galaxies for which photometric redshifts are 
available in the literature (Coe et al. 2006) , and present the 
mean scalelengths for each magnitude and redshift bin in 
Table 1 (see also Sec. 3.1). 
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Figure 3. Normalized distribution of scalelengths three absolute 
magnitude and redshift bins. The brightest absolute magnitude 
bin corresponding to the orange and red data points in Fig. 2 are 
shown here. The SDSS data are illustrated by the black histogram 
(10040 galaxies) and HST data in red and blue histograms with 
125 and 26 galaxies, respectively. Similar analysis performed on 
the photometric redshift galaxies yields the mean scalelengths 
presented in Table 1. 



3 OBSERVATIONAL EFFECTS 
3.1 Biases 

At the highest redshifts, the fight observed in the z— band 
will have been emitted as ultraviolet photons which are pro- 
duced by hotter stars. Ideally, a commensurable compari- 
son of all scalelengths would have to account for wavelength 
shift due to a galaxy's redshift. Serious problems arise how- 
ever, since at each redshift, we would be required to have 
detailed knowledge about dust properties, stellar popula- 
tions, and galaxy evolution models for each morphological 
type. Such models are still highly uncertain at low redshift. 



and unavailable for high redshift galaxies. Studies of galax- 
ies below redshift 1.0 have found that for a typical 10^^ M0 
galaxy, the salelength is ~ 20% larger in the B— band as 
compared to the A'— band (Peletier et al. 1994; Giovanelli 
& Haynes 2002). Thereof, wavelength dependencies cannot 
account for more than about 20% of scalelength differences 
found here. Moreover, increased dust content makes a host 
galaxy more opaque which leads to an artificially larger bluer 
band scalelength when compared to red and/or infrared im- 
ages (Valentijn 1990). For a subsample of our 686 galaxies, 
for which reliable HST- ACS images with the F435W fiher 
(B— band) are available, we have derived the scalelengths 
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Figure 4. Scalelength versus photometric redshift for the 
HUDF (coloured points) and GOODS galaxies (open circles). All 
GOODS galaxies above redshift 2.0 have absolute magnitudes be- 
low -22. 



four smaller than their brightest counterparts. In the red- 
shift range 2.0-5.8 galaxies in all these luminosity ranges 
have scalelengths that are significantly smaller than their 
local counterparts. In the intermediate redshift range, the 
Universe is populated by galaxies with intermediate scale- 
length values (see Table 1). The photometric redshift sample 
further shows that scalelengths remain redshift-invariant for 
systems fainter than absolute magnitude > — 18 and below 
redshift 1.5. 

High redshift galaxy population may differ significantly 
from galaxies in the local Universe. The effects of stellar 
population and dust content on the measured scalelengths 
are of the order of 20%. The effect of spatial resolution and 
sampling accounts for an additional 15%. Adding all these 
effects quadraticly, the resulting scalelengths at high redshift 
are ~ 35% lower than their intrinsic values. 



at both wavelength ranges. We found that at redshifts be- 
tween 2.0 and 5.8 the z— band scalelengths are smaller than 
those measured in B— band by ~ 20%. This indicates that 
the combined effects of stellar populations and dust add up 
quadraticly. 

At the redshift range between 2.0 and 5.8 the mean disc 
scaelength is 0.6 ± 0.2 kpc, corresponding to 2-3 HST-ACS 
pixels. Although the scalelength is not an actual galaxy size, 
it is imperative to quantify the effect of smearing on our de- 
rived scalelengths at high redshift. To ensure that the mea- 
sured scalelengths do not contain any systematic redshift 
related bias due to differences in resolution for each redshift 
slice, we have artificially redshifted a sub-sample of 23 low- 
redshift galaxies, and re-derived their scalelengths at each 
epoch. All test galaxies were located at redshifts below 0.5, 
and we simulated their appearance at redshifts 2, 3, 4, and 5 
using the redshift simulation code developed by Epinat et al. 
(2010) that simulates the reduction of resolution and rescal- 
ing of the point spread function of each image at the desired 
redshift. We found that, at the highest redshift, poorly re- 
solved galaxies exhibit artificially lower scalelengths by at 
most 15% of the initially observed low-redshift values, and 
Sersic indices are artificially lowered by ~ 20%, whereby the 
presence of an exponential component is re-confirmed. 

We further test the effect of cosmic variance, and find 
the same redshift-scalelength relation for each individual 
survey used here. We test the effect of Malmquist bias and 
include galaxies for which photometric redshifts are avail- 
able. These galaxies are typically fainter, however their red- 
shift measurements are not as reliable as we would require 
for our analysis, nevertheless, they still give us useful clues 
about a possible Malmquist bias effect. Deriving scalelengths 
and applying our selection criteria to a larger sample of 1472 
galaxies with photometric redshifts from the HUDF (Coe et 
al. 2006), we find similar trends based on the distributions of 
the photometric redshift galaxies (see Table 1). We further 
note that the scalelength decrease in the GOODS galaxies 
(Bundy et al. 2009) agrees with that found for the HUDF 
sample (see Fig. 4), however, as we opt to use a homoge- 
nous set of photometric redshift data set and as the HUDF 
contains more galaxies at redshift > 2.0 (for which we have 
reliable scalelengths), the values in Table 1 are based on the 
HUDF galaxies only. 

In the local Universe, galaxies with absolute magni- 
tude between —18 and —20 exhibit scalelengths a factor 



3.2 Surface Brightness Limitation 

Due to the degeneracy between surface brightness selec- 
tion effect and direct luminosity and size measurements, 
the HST images analyzed here would not allow direct de- 
tection of galaxies with half-light radii above ~ 25 pixels 
at 2— band AB magnitude fainter than 24, and above ~ 15 
pixels at z— band magnitudes fainter than 26. As found by, 
e.g., Bouwens et al. (2006), significant loss of the large low 
surface brightness galaxies are expected to affect the data 
at z— band AB magnitudes fainter than ~ 27 and axial ratio 
analysis of disc galaxies out to redshift 4 has revealed that 
only 25% of the spirals larger than 10 pixels are lost below 
the HST detection limit (Elmegreen et al. 2005b). These lim- 
its are significantly fainter than the z— band AB magnitude 
limit of our spectroscopic and photometric redshift sample, 
~ 24 and ~ 26, respectively. 



4 DISCUSSION 

Our results agree with previous works postulating that the 
first generation galaxies are expected to be gas-rich, and 
when they undergo mergers, they leave behind substantially 
rotating remnants and thereby galaxy discs are formed early 
in the Universe (Wuyts et al. 2010). The discs then continue 
to grow to transform into the Hubble sequence observed in 
the local Universe. Comparing our derived scalelengths with 
those measured by means of direct size measurements, we 
find that direct size measurements typically deliver over- 
estimated values due to their strong dependancy on the ob- 
served luminosity in the sense that they do not pick out the 
larger lower surface brightness galaxies (e.g., Bouwens et al. 
2004). Further advantage of deriving scalelengths as opposed 
to direct size measurements for high redshift galaxies can be 
seen in the comparison with predictions from cosmological 
galaxy formation simulations with the computational reso- 
lution that allows the calculation of disc scalelengths (e.g., 
Ceverino et al. 2010; Dutton 2009). On the observational 
side, to date, scalelength measurements have been usually 
based on inhomogeneous data sets from the local Universe. 
It is only two years ago that scalelengths were derived for 
a comprehensive sample of 30000 SDSS galaxies out to red- 
shift 0.3 (Fathi et al. 2010), while studies beyond redshift 
~ 1.5 were limited to a handful of objects (van der Kruit 
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Scalelength (kiloparsoc) 
< redshift < 0.3 


Scalelength (kiloparsoc) 
0.3 < redshift < 2.0 


Scalelength (kiloparsec) 
2.0 < redshift < 5.8 


(-24,-22) 


4.8 ± 1.5 


10040 galaxies 




2.0 ±0.6 


146 galaxies 




0.7 ± 0.3 


56 galaxies 




(-22, -20) 


2.7 ± 1.2 


17802 galaxies 




1.3 ±0.5 


514 galaxies 




0.7 ± 0.2 


110 galaxies 




(-20,-18) 


1.2 ± 0.4 


1886 galaxies 




0.8 ±0.3 


281 galaxies 




0.6 ±0.2 


60 galaxies 





Table 1. We test the validity of our results based on spectroscopic redshifts by expanding the sample to include galaxies for 
which photometric redshifts are available. All values presented here are based on Gaussian fits to the V/Vmax volume— corrected 
scalelength distribution of the SDSS sample combined with the 1472 photometric redshift galaxies in the HUDF. In each cell, the 
number of galaxies are written at the right. The errors are the standard deviations of the Gaussian fits (c.f. Fig. 3). The weighted 
means remain within the uncertainties when using both HUDF and GOODS galaxies (c.f.. Fig. 4). 



& Freeman 2011). Here we have compared our results from 
the SDSS sample with scalelengths derived from a sample of 
686 galaxies with spectroscopic redshifts 5.8. While we use 
1472 photometric redshift galaxies to test against Malmquist 
bias, our results are based on solely spectroscopic redshift 
measurements, whereas all previous works dealing with size 
evolution of disc galaxies, have been based on photometric 
redshifts or drop-out techniques. 

Our results confirm the presence of a dominating ex- 
ponential components in galaxies out to redshift 5.8 (c.f., 
Ravindranath et al. 2006), and we find that at redshifts 
between 2.0 and 5.8 (corresponding to approximately 10- 
12.5 Gyr ago), bright disc galaxies had scalelengths up to 
a factor 8 smaller than their local counterparts. A similar 
trend for galaxies fainter than absolute magnitude —22 can- 
not be confirmed. The observed decrease out to redshift 5.8 
is significantly greater than the value predicted by the cos- 
mological picture in which baryonic disc scalelength scales 
with the virial radius of the dark matter halo. At the highest 
redshifts, direct size measurement methods suggest a rather 
slow disc size evolution, whereas our results show a disc size 
growth speedier than that predicted by the H~^^^ model 
(e.g., Bouwens et al. 2004), and slower than that predicted 
by the model (e.g., Ferguson et al. 2004). Our findings 
suggest that the bulk of galaxy growth occurred in the first 
1-3 Gyr after Big Bang, i.e., right after the formation of the 
first galaxies, and that a combination of both models should 
be considered in cosmological galaxy formation models. 

The observational picture presented here offers a 
uniquely quantified and newly examined scenario for mod- 
ellers of galaxy formation and evolution within a cosmologi- 
cal paradigm. Applying a similar analysis to wide area high 
redshift surveys such as CANDELS (Grogin et al. 2011) will 
help building larger sample and improve measurement of 
galaxy structural parameters at high redshifts. 
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